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1.  Introduction 

The  vibrations  of. a  string  of  density  p(x),  length  L  taut  by  a  unit 
tension,  are  governed  by  the  equation 

u"  +  w2pu  =  0,  xe(0,L).  (1) 

2 

It  is  well  known  that  if  peC  (0,L),  this  equation  can  be  transformed  into 
the  canonical  Sturm-Liouville  equation 

2 

y"  +  (w  -q)y  =  0. 

However,  even  in  those  cases  in  which  p  is  sufficiently  smooth,  there  are 
certain  advantages  in  dealing  with  the  problem  in  its  original  form  (1), 
particularly  if  one  is  interested  in  other  vibrating  problems. 

In  the  present  paper,  I  shall  present  a  new  approach  to  the  solution 
of  the  inverse  problem  for  the  vibrating  string,  i.e.  to  the  reconstruction 

CO  CO 

of  the  density  p(x)  given  its  length  L  and  two  spectra  and  • 

In  order  to  simplify  the  presentation,  I  shall  assume  that  the  spectrum 
{>.  }"  is  associated  with  (1)  and  the  boundary  conditions 

u(0),  =  u(L)  =  0  ,  (2) 

00 

i.e.  the  fixed/fixed  configuration,  while  is  associated  with  the  free/ 

fixed  configuration: 

u'(0)  =  u(L)  =  0.  (3) 


3 


2 

u"  +  u>  pu  =  0  , 
u(£)  =  u(L)  =  0  , 


00 

and  {yn(£)}j  the  spectrum  of 


2 

u"  +  w  pu  =  0  ,  ^ 

u'(£)  =  u(L)  =  0  ,  ( 


where  0  <_  i  <  L.  If  p(x)  is  continuous  and  bounded  away  from  zero  for 
xe(0,L),  then 


pW  = 


w 

IT 


<L-*>  *i  W  n-1  »».l  «“„<*> 


Theorem  2:  Let  {pn>~  and  L  be  such  that  the  solution  p(x)  of  the 

inverse  problem  for  the  vibrating  string  exists  and  is  continuous  and  bounded 
away  from  zero.  Then,  this  solution  is  given  by 


P(x)  = 


(L-x)2  P12(x) 


*n4(x) 


n=l 


P_.,2(x)  Pn2(x) 


n+1 


00  00 

where  and  {un(x)}j  are  the  solutions  of  the  initial  value  problems 


dx 


»„2w  ij, 

L-x  »  — 


"["[  Cl-Xn2(x)/Xk2(x)) 


kj*n 


-t; 


du  -  .  (1-Un2(x)/Xk2(x)) 

-dir  -  Mn4(*)p(x)(L-x)  r - 


kj*n 


U-Pn2(x)/uk2(x)) 


with 


=  *n2  * 


2  2 
“n  <°>  '  »„  ' 


3.  Proofs 


Proof  of  Theorem  1.  The  given  data  are  such  as  to  insure  the 
existence  and  uniqueness  of  the  solution  of  the  inverse  problem.  Following 
Krein  [1,2],  we  can  construct  this  solution  by  considering  a  converging 
sequence  of  approximations  to  the  integrated  mass,  viz. 


rx 


m(x)  = 


p(x')dx' 


.(N), 


The  generic  Nth  order  approximation  to  m(x),  say  mv  (x) ,  is  a  piece-wise 
constant  function: 


m 


(N) 


X 

(x)  =  l  m.  for  xe[x. ,  x.  .),  i=l,2,...,N 

3=0  J  lii 


(4) 


In  the  above  formula  m  =0,  x  =0  and  x..  ,=L.  The  other  values  of  m.  and 

oo  N+l  1 


x. ,  or  rather 
i 


•  v, 


are  found  by  writing  the  rational  fraction  -Uo(u>  )/U0'(u>  ) 

N 

Uo(u)2)  =  -L  |  |  (1-  o>2/Xn2)  , 

n=l 


and 


U0'(„2) 


N 

TT ci-  u2/v 2)  , 

n=l  n 


as  a  Stieltjes  continued  fraction,  viz. 


V“2) 

uo'(«2) 


l  + 
O 


It  is  possible  to  write  down  explicit  formulas  for 

2  2 

To  that  effect,  we  first  divide  -UQ(a>  )  by  Uo'(u  ),  viz. 


-uoo2)  =  *0u0'(»2)  '  ui(a,2)  ’ 


2 

and  then  LMfu)  )  by  the  remainder  of  the  previous  division, 
-UjCoj2): 

U  '  (oj2)  =  -mjU^jCu2)  +  Uj '  (u>2)  . 


where 


l  and  m, 
o  1 


namely 


■_1*  ^  *  k  '1. 


,d  Vs. 


6 


Clearly 


N 


(9) 


The  polynomial  U^(o)  )  being  of  degree  (N-l)  can  be  written  as  follows: 


N-l 

Uj(to2)  =  -(L-Ao) 

n=l 


(10) 


where  X^  is  the  nth  eigenfrequency  of  a  string  with  point  masses 
N  N 

(m. }  located  at  (x.}»  whose  left  end  at  x=x,h£  and  right  end  at 
1  /  1  /  1  o 

2 

x=L  are  held  fixed.  With  this  expression  for  )  we  can  carry  out 

the  division  in  (8)  and  get 


ml  =  (L-Xj) 


N-l 


n=l 

N 


n=l 


Now  consider  the  ratio  m,/£  :  by  definition 

l  o 


£ 

o 


(x)dx 


x 


1 


(11) 


(N) 

where  7 (x)  is  the  N-th  approximation  to  the  density.  If  p  is  bounded 
away  from  zero,  then  as  Nx»,  x^-*-0  and  the  above  ratio  tends  to  p(0). 


Therefore 


7 


P(0) 


Aim 

N-*=° 


1 

(L-Xj) L  yj2 


N-l 


n=l 


2  2 
Vl 


As  N-**>,  Xj+0  and  consequently  ;  making  use  of  the  asymptotic  form  of 

the  eigenvalues,  we  can  show  that  the  resulting  infinite  product  is  con¬ 
vergent  and  consequently  that 


p(0) 


.2  2 
L  U1 


Vl  Pn 


n=l 


(12) 


By  repeating  the  analysis  for  that  portion  of  the  string  occupying 
the  segment  (£,L)  we  could  deduce  that 


pW 


_ l _ 

(L-A)2  yj2 (A) 

n=l 


(*) 


2 

vn+l 


(A)yn2(£) 


(13) 


Proof  of  Theorem  2. 
(1)  such  that 


Let  U(x,u>  )  denote  a  fundamental  solution  of 


U(L,u>2)  =  0  , 


U'(L,u)2)  =  1  . 


8 


2  2  2 
Since  U(x,u>  )  and  U'(x,to  )  are  entire  functions  of  to  of  order  1/2,  they 

can  be  written  as  follows: 


U(x,o>2)  =  -(L-x)  |  |  (1-J/X  2(x))  , 

n=l 


U  *  (x,'u)  )  = 


I  |  (l-to2/y  2(x))  . 


The  coefficients  in  front  of  the  infinite  products  have  been  determined  by 

2  2  2 
setting  a)  equal  to  zero  in  (1).  The  zeros  of  U(x,a>  )  and  U'(x,w  )  are 

the  eigenvalues  of  that  portion  of  the  string  occupying  the  interval 

(x,L)  and  vibrating  in  the  fixed/fixed  and  free/fixed  configurations. 

The  product  representations  (14)  must  (i)  be  compatible  and  (ii) 

satisfy  the  original  equation  (1).  As  a  result,  we  must  have 


(l-0)2Ak2(x)) 


-  (L-x)  l 


n=l  X  (x) 
n  v 


(1-  - ) 

Xk  (x) 


(1-  -f -  )  , 

v« 


»  2  dp  2 

y  _ D_  a _ “ _ ) 

n=l  u  4(x)  dx  w/(x) 

n  k/n  K 


=  to2p(x)(L-x)  (1-  — | - ) 


1  V<x> 


9 


2 

By  equating  the  coefficients  of  equal  powers  of  w  ,  the  above  identities 
give  rise  to  an  infinite  system  of  linear  equations  for  the  unknowns 

2  oo  2  00 

(dX  /dx)  and  {dp  /dx)  .  The  solution  of  this  infinite  system  can 
n  1  n  1 

be  obtained  via  the  theory  of  infinite  determinants.  However,  there  is 
a  simpler  way  of  arriving  at  the  solution.  To  that  effect  we  return  to 
(14)  and  write 

U(x,Xn2(x)  =  0  ,  n=l,  2,  ... 


Differentiating  with  respect  to  x,  we  get: 


9U(X,h)  ) 
9x 


dX 


2  2 
«  =  A/(x) 


dx 


9U(x,ta  ) 
9u2 


=  0 


0)  = 


*>) 


and  making  use  of  the  product  representations  (14),  we  deduce  that 


Similarly 


dX_ 


dx 


X  2(x) 
n 

L-x 


X  2(x) 

(l-  A 


k=I  « 


) 


X  2(x) 

(i-  A — ) 


(15) 


du 


dx 


V*n  (x)  P  (x)  (L-x) 


yn2(x) 

Cl-  A — ) 


U  *k  w 


k^n 


Vn2(x) 

(i-  A — ) 

"k  w 


(16) 


10 


Or,  substituting  p(x)  by  the  expression  given  in  (13), 


dx 


1 

L-x 


V* !  2  (x) 


A  4(x) 
m  v 


m=l 


|  (1-Pn2(x)/Xk2(x)) 

J^l _ 

00 

|  |  U-Pn200/uk2(x)) 

k^n 


(17) 


We  shall  establish  next  that  the  initial  value  problem  consisting  of 
(15),  (17)  subject  to  the  initial  conditions 

(n=l ,  2,  ...)  (18) 

has  a  unique  solution.  For  convenience,  we  can  write  this  initial  value 
problem  as 


2  2 
A  z(0)  =  X 
n  K  J  n 


Pn2(°)  -  F„2 


^=h(v)  , 

v(0)  =  v  , 

where  the  stretched  coordinate  ?  is  defined  thus 

C  -  In  (1  -  £  )  , 

v  stands  for  the  sequence  obtained  by  interlacing  the  eigenvalues 

(p2)"  and  {X2}?  ,  namely 
n  1  n  l 


12 


Since 


.  2  2  2.,(l 

Kn  “  n  *  /CJ 


pW 


Pn2  -  p^dt)2 

Jq 


for  n-**> 


it  is  sufficient  to  show  that 


F  = 

n 


e 


i  -inn 
n 


1/k 


=  0(1) 


and 


1  ,  U 

r  „  ---inn  *n 
G  =  e  n  — r- 

n  2 

*1 


.  2  2 

”^k  -1/k 

0  +  -ky- kx-  )  e  1/K  -  0(1) 
"k  -wn 


k^n 


As  a  matter  of  fact,  we  shall  only  show  that  (21a)  holds;  the  same 
arguments  can  be  used  to  prove  (21b) .  Since 


for  all  Mn  , 


we  can  bound  F  above  and  below: 
n 


13 


1  !  A  2-y  2 

exp  [  -  -  inn  +  J  (  k  -  )]  1  F 

k/n  y.  -A 

K  n 


5  exp  [  i  -  inn  +  £  (  J- 


k^n 


,  2  2 
k  ‘"k 

k  n 


For  large  values  of  n,  say  n>N,  where  N  denotes  the  subscript  above  which 
2  2 

A  and  u  take  on  their  asymptotic  values,  we  can  rewrite  the  above  in- 
n  n 

equalities  as: 


„  wi2  2  „  ,  ,  2  2 

n  ,  NX,  -y,  n-1  X,  -y, 

exp  (Cl  f  -  «n„)  .  £  *  1  ~T^T 

1  IX  -y.  N+l  X  -y. 

n  k  n  *k 


2  2  2 

T  i  vV 

I  (r- 


,  2  2 


n+l 


2  ,  2 
y,  -X 
k  n 


>  *  l 


4n2+l 


y,  -X 
k  n 


ii  ^ 

£  exp  [ (  l  7  -  *nn)  +  l  ?  2 


«.  .2  2  ,  ,  2  2 

N  X.  -y  n-1  X  -y 
K  k  k  k 


1  X  -A, 
n  k 


2  2 
N+l  X  -A. 
n  k 


4n  1 

*  I  Cr 


,2  2  %  2  2 

X,  00  .  X,  -u. 

k  k  )  +  I  (f  - 


-  ,2,2'  %  vk  ,2,2 

n+l  X,  -X  .  2x,  X.  -X 

k  n  4n  +1  k  n 


)] 


Now 


N  - 

l  «  0(n-Z) 


14 


l  =  j  Ann  +  o(l) 
N+l 


1 

l  =  -  -  £nn  +  o(l) 
n+l  * 


l  =  o(l)  , 

4n2+l 


where  we  have  used  the  sigma  sign  as  a  shorthand  notation  for  the  various 
sums  entering  in  the  exponents.  Therefore, 


v 

F  ->  e  as  n-*» 


where  y  is  Euler's  constant. 

The  next  step  consists  in  introducing  the  sequence  space  S  where 


S  = 


00 


IMI  -  l 

1 


(22) 


as  well  as  the  supspace  S*  C  5  of  sequences  which  (i)  are  positive  (ii)  are 
increasing  and  (iii)  satisfy  the  asymptotic  relation 


2  2 

n  tt 

v  ~  - =•  as  n+® 

n  4I2 


(23) 


where  1  is  a  constant.  Clearly,  sequences  obtained  by  interlacing  two 
spectra  lie  in  S*.  Also,  h  can  be  thought  of  as  a  map  of  S*  in  S. 


15 


The  spaces  S  and  S*  are  convex:  indeed,  if  o,neS*,  then 
ao+(l-a)neS*  where  0 <a <1 .  More  importantly,  the  functions  hn(ao  + 
(l-a)n)  ,  n=l ,  2,  ...  are  analytic  functions  of  a  for  ae[0,l].  Con¬ 
sequently,  for  every  n  there  exists  a  value  afi  such  that 


hn(o)  -  hn(n) 


dh 
_ n 

da 


-j=a 

n 


or  equivalently 


3h 

-  \Cn)  -  l 

m  d  a 

m 


^°2m  2m^ 


3h 

+  l  —~ 

m  3y  1 
m 


C°2m-rT12m-l) 


We  should  remark  at  this  stage  that  on  account  of  the  asymptotic  form  of 
h^,  the  an's  become  independent  of  n,  or  more  precisely 


2n 


2n-l 


as  n-*» 


2 

Referring  to  (191,  we  can  deduce  explicit  formulas  for  3h  / 3A-  and 

n  in 

2 

3h  /3u  .  As  a  result 

n  m 


16 


h2n(0)-h2n(n)  = 


-X 


(2n)  ,(2n) 


n 


l 

mj*n 


°2m  '  n2m 


x(2„)2  (2r,)2  C2H)2, 

m  m  n 


+  f^2n)  [ 
n  1 


(2n)‘ 


-  I 


k  (2n)2_A(2n)2 
yk  n 


kjta  A(2n)\x(2n) 
Ak  n 


2^ (a2n"n2n) 


♦  x<2")2  f„ 

n 


(2n) 


°2m-l  "  n2m-l 


m  (2n)2  (2n)2_x(2n)2 

vm  lum  n 


(24a) 


where 

a2u  °2m  *  (1-“2n)n2«  ' 

°2„  °2-l  4  (*-“2n,n2m-l  1 

and 

fn2"’  '  h2n  l“2n°  *  ' 


x(2n)2  , 
m 

(2n) 2 
% 


Similarly 

h2n-l(<”-h2n-l 


(n)  =  n 


(2n-l) 

n 


2  (2n-l) 
ftn 


l  °2.  ~  n2m _ 

■  x(2"-i>2c)l<2n-l)2  (2n-i)2 

m  v  m  n 


_ii(2n-l)2 


(2n-D 

6n 


l 

m/n 


°2m-l  ~  n2m-l _ 

(2n-D2  (2n-l)2  (2„-l)2 

v  m  n 


I 


+  2g^n_1)  ][  {  2m  n2m  _  2m-l  ~  n2m-l  ^ 
n  m  x(2n-l)2  (2n-l) 2 


(24b) 


+  (2n-l)  .  2  r  _ 1 _  r.  1 

8n  (2n-l)2  •  k  A(2n-1)2  (2n-l)2  k/n  (2n-l)2  (2n-l)21 (°2n-1  ' 


Because  of  the  asymptotic  form  of  f  and  g  ,  we  can  assume  that  h  (o)-h  (n) 

n  n  n  n 

is  positive  for  n>2N,  by  interchanging  if  need  be  a  and  n.  Dividing  (24a) 

4  4 

by  (2n)  and  (24b)  by  (2n-l)  and  summing  over  n,  we  can  see  that 


v  hn(c)_hn(n) 

V  4 
1  n 


or,  in  view  of  the  previous  remark 


<_  const.  ||o-n||  , 


|h(cr)  -h(n) 


,2N  IV«)-hnO) 

1  „4 


<_  const.  1 1 a-rj | 


2N  h  (a)  -h  (n) 
y  n  n 

f  4 

L  n 


In  order  to  derive  the  above  inequality,  we  have  made  use  of  the  fact  that 
the  dependence  of  the  various  quantities  on  can  be  ignored  for  large  n's 
and  that  the  sums 


y  _ I _  y _ 1 _  y  1  y  _  1 

4.  2  ,  2  '  1  .  2  2  •  L  2.2*  4  2  2 

n?m  >  -X  n  X  -u  n  p„  -X  npm  p  -p 

n  m  n  m  n  m  n  m 


..  . -  -r 
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-2 

are  of  0(m  ).  Noting  that  (24a)  and  (24b)  imply  that 


lhn(o)-hn(n) |  <  const.  ||o-n|| 


(26) 


we  therefore  conclude  that 


I  |h(a)  -  h(n)  1 1  <_  const.  ||a-n|!  • 


(27) 


This  Lipschitz  condition  guarantees  the  existence  of  a  unique  solution 
to  the  initial  value  probii^. 

4.  The  Quasi -homogeneous  String 

This  paragraph  examines  what  happens  to  formula  (13)  for  the  case 
of  a  quasi -homogeneous  string,  i.e.  for  a  string  whose  density  is: 


p  =  p  [1  +  er  (x)  ] 


(28) 


where  e<<l.  Making  use  of  perturbation  theory,  one  can  easily  show  that 


and 


2  2 

»  -Jn~-  -2  f1  ♦  E(-r0  +  J  r2n}  +  0(e2)1  * 


n 


p(L-£) 


.2  2 


"n2«  ■  l1  -  *  1  r2„-l’  *  0(e2)1  - 

p(L-£) 


(28a) 


(28b) 


where 


=  L-*  J 


r(x)dx  , 


(29) 
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and 


r 


n 


2 

L -2 


t  rM 


2 


COS  H7T 


x-2 

L-2 


dx 


(30) 


Therefore 


TTTT 2  11  *  c(r2n  *  7  T2n.l  *  7  r2n-l’  *  0('2)1 
(n  -h) 


and  consequently 


00  00 


4 

X  (2) 
n 

4 

n 

n= 

u2  ,(2)u  2 (2) 
n+1  JVn  v  J  i 

i= 

,  2  i^2 

(n  -h) 

[1  *  E  j  (r2n 
n=l 


r  +r 
2n+l  2n-l 

2 


■)*  0(e2) ] 


2  r  00 

V  (1  ♦  e(-  4+  X  O  +  °Ce2)>  • 

*  n=l 


Hence 


Cl-t)V W  l=[  n2(»)  '  ,  .  e(r  *  ^  *  OC.2) 


\>) 


r  c»  - 

1  +  e(*T4[rn)  +  °CE  ) 


=  P  [1  +  e  [  rn  +  0(e)] 
o 


(31) 


Thus,  for  the  simple  case  of  the  quasi-homogeneous  string,  formula  (13) 
reduces  to  an  end-point  evaluation  of  the  Fourier  cosine  series  for  r(x) . 
Indeed,  from  the  representation 


20 


r(x)  =  r  + 
o 


cos  rnr 


x-fc 

i-a 


it  is  obvious  that 

oo 

p(£)  =  p  [1  +  e  l  rn]  . 

o 

5.  Concluding  Remarks 

The  derivation  of  p(0)  via  the  computation  of  JIq  and  m^  can  be 
thought  of  as  involving  a  process  in  which  the  frequency  w  tends  to  infinity. 

Had  we  started  with  a  time  dependent  formulation  of  the  problem,  this 
process  would  have  been  associated  with  letting  the  time  t  tend  to  zero. 

In  that  sense,  the  method  is  reminiscent  of  that  of  Kac  [3],  Deift  5 
Trubowitz  [4]  and  Trubowitz  [5] .  A  similar  approach  can  also  be  used  for 
solving  the  canonical  Sturm-Liouville  problem:  in  this  case  q(x)  is  expressed 
via  a  trace  formula  such  as  the  ones  given  in  [6,7,8].  However,  the  results 
are  not  as  explicit  as  in  the  present  case. 

Finally,  we  have  restricted  our  analysis  to  the  boundary  conditions 
(2)  and  (3).  If  these  boundary  conditions  are  replaced  by  the  following 
ones: 

u(0)  =  u(L)  cosy  +  u'(L)siny  =  0  (1') 

and 

u(0)  cosg  -  u’(0)  sing  =  u(L)  cosy  +  u'(L)  siny  =  0  (2') 


where  0  <  0  ir/2  and  0  <^y  <  ir/2,  then  the  expressions  for  and  nij 
change  slightly  (see  [9])  and  as  a  result  formula  (13)  becomes: 


p(x)  = 


sin2B  (L-x)2  P12(x) 


*n4(x) 


n=l 


il^x)pn2(x) 


(13’) 
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